Quantum atom optics with fermions from molecular dissociation 
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We study a fermionic atom optics counterpart of parametric down-conversion with photons. This 
can be realized through dissociation of a Bose-Einstein condensate of molecular dimers consisting of 
fermionic atoms. We present a theoretical model describing the quantum dynamics of dissociation 
and find analytic solutions for mode occupancies and atomic pair correlations, valid in the short 
time limit. The solutions are used to identify upper bounds for the correlation functions, which are 
applicable to any fermionic system and correspond to ideal particle number-difference squeezing. 
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Advances in the experimental control of degenerate 
quantum gases of neutral atoms have recently reached 
the stage where atomic correlations and quantum statis- 
tics can be directly accessed via the measurement of atom 
shot noise and atom counting [lj, |2|, Li Lit l|- Given 
the similarities with the pioneering photon correlation 
measurements of Hanbury Brown and Twiss [|J and the 
intriguing parallels with modern quantum optics, these 
experiments represent a remarkable step forward in ad- 
vancing quantum atom optics. Earlier and related exper- 
iments were performed either in a cold but not degener- 
ate atomic beam |j| or else were examples of indirect 
measurements of higher-order correlations and number 
squeezing |8l Ifl Ulfllllj . 

In quantum optics with photons the most successful 
applications have been achieved using squeezed light and 
entangled photon pairs from parametric down conversion 
[l2l Hil Il4 Hif . A matter- wave or atom optics counter- 
part of down conversion can be realized through disso- 
ciation of a Bose-Einstein condensate (BEC) of molec- 
ular dimers. The most direct analogy with quantum 
optics corresponds to the case of bosonic statistics of 
constituent atom pairs 0, 0, 0, 0], as realized in 
dissociation experiments with 23 Na2 and 87 Rb2 |2(ll2l|. 
In contrast to this, the recent correlation measurements 
at JILA 0], using dissociation of 40 K2 molecules near a 
magnetic Feshbach resonance, have an intriguing twist in 
that the constituent atoms obey fermionic statistics. In 
this case, analogies with photonic down conversion and 
the implications for possible future applications are not 
so immediate as in the case of bosons. 

In this paper we study the quantum dynamics of dis- 
sociation of diatomic molecules into fermionic atoms and 
analyze the resulting atom correlations in momentum 
space. At low densities and small s-wave interactions the 
momentum correlations are reflected in the density corre- 
lations after spatial expansion of the cloud. Accordingly, 
the results obtained here are related to the spatial corre- 
lation measurements performed at JILA pj. Together 
with the recent studies of fermionic four-wave mixing 
l'2/fl and association of fermionic atoms into molecules 



[23I Eif . the present work (see also [25|) expands the 
paradigm of fermionic quantum atom optics. Possible ap- 
plications are in precision measurements and fundamen- 
tal tests of quantum mechanics, similar to those proposed 
recently for dissociation into bosonic atoms [13, LUj • 

Due to the limitations of the undepleted molecular 
field approximation employed here, the obtained results 
are only applicable to short dissociation times. How- 
ever, the advantage is the analytic transparency of the 
results, which provide useful insights at the conceptual 
level. In particular, we show rigorously that the obtained 
pair correlations represent their generic upper bounds 
and are applicable to any fermion system. This gives 
a useful reference for further (numerical) studies of this 
and related systems with less restricted approximations. 
In addition, we point out that the notions of maximum 
correlation and squeezing of atom number-difference fluc- 
tuations have to take into account the fact that the shot 
noise level for fermions is fundamentally different to what 
one usually encounters in quantum optics with bosons. 

We start the analysis by considering an effective field 
theory Hamiltonian for the coupled atomic-molecular 
system given by H = i? - ih\ J dx(*J* T *| - *j\&{\& ) 
|2^. Here, Ho stands for the usual kinetic energy term 
plus the trapping potential, ^(x, i) is the bosonic field 
operator for molecules of mass mo, while ^jQ^x, t) are 
fermionic operators for atoms (with masses nj|jn = 
mo/2 = m) in two different spin states, a =t,4- The 
atom-molecule coupling is described by and we have 
omitted intra- and inter-species s-wave scattering inter- 
action terms, which is justified at low particle densities. 

Considering a uniform system in a cubic box of side L, 
we employ an undepleted molecular field approximation 
valid in the short time limit. The molecular field is de- 
scribed by a coherent state and we absorb its mean-field 
amplitude ^0 into an effective coupling g — Xa/^O' wnere 
no = \^o\ 2 is the density. Assuming periodic boundary 
conditions and expanding the atomic fields in a plane- 
wave basis, in terms of single-mode operators, we obtain 
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the following effective Hamiltonian, in a rotating frame: 
H = J2 k a ^k^ka - i%^ k (c kT c_ ki - cl ki c kT ). (1) 

Here, c ko . (ck<r) are fermionic creation (annihilation) op- 
erators, fiua = c k(T Cko- is the particle number operator, k 
is the momentum [fej = 2Ttrii/L, rii — 0, ±1,±2, 
z = x, y, 2] and A k = A + hk 2 /(2m). The detuning 2A 
corresponds to the overall energy mismatch 2hA between 
the free two-atom state in the dissociation threshold and 
the bound molecular state. 

The model system in mind corresponds to a pure 
molecular BEC on the stable side of a Feshbach reso- 
nance, with no residual atoms present. This is followed 
by a rapid switching on of the coupling \ (e.g., via an rf 
transition or a rapid crossing through the resonance to 
the atomic side) and a simultaneous switching off of the 
trapping potential. From this stage onward, the atomic 
filed evolves in free space, with a vacuum initial state and 
negative detuning A. Since the excess of energy 2fi|A| is 
released into the kinetic energy of dissociated atom pairs 
in the two spin states, 2fr|A| -> fr 2 (|k T | 2 + |kj 2 )/(2m), 
we expect - from momentum conservation - strong cor- 
relation between the atoms with opposite spins and mo- 
menta, kf = — kj. In fact, the interaction term in the 
Hamiltonian Q is the prototype interaction to produce 
- in the lowest order perturbation theory - an entangled 
spin singlet state. It is also a fermionic analog of the 
squeezing Hamiltonian in quantum optics [l5| . 

Introducing a dimensionless time t = t/to, length I = 
L/do, detuning S = Aig, and momentum q = kcfo, where 
to = l/ff is the time scale and do = y/hto/ (2m) is the 
length scale, we can put the system into a dimensionless 
form, with the Heisenberg equations of motion 



dc^/dr = -i<S q c qT - ct qi , 



(2a) 
(2b) 



where <5 q = q 2 + S (q = |q|). Solutions to Eqs. (J2J) are: 

c qT (r) = A q (r)c qT (0) - S q (r)cl qi (0), (3a) 
5l qi (r) = S q (r)c qT (0) + A q (r)c! q; (0). (3b) 

Here, A q (r) = cos (g q r) - i<5 q sin (p q r) /g q , B q (r) = 
sin( 5q r)/ 5q , ,g q = (1 + tf 2 ) 1 ^ an d |A q | 2 + B\ = 1. 
Using these solutions with vacuum initial conditions, we 
find that the only nonzero second-order moments are the 
mode occupancies and the pairing fields: 

n q (r) = (n qCT (r)) = B\{t) = sin 2 ( 3q r) / 5 2 , (4) 
m q (r) = (c^^c^t)} = A q (r)B q (r), (5) 




which, in addition, are related by 

l"t q (r)| 2 = [l-n q (r)]n q (r). 



(6) 



FIG. 1: (a) Average occupancies of one of the resonant modes 
lo = ypf and two sidebands q t (full lines) as a function of 
time, for S — — 16. The dashed line is the solution for the go- 
mode in the case of bosons, which grows exponentially. The 
inset shows the slice through the origin of the fermionic 3D 
momentum distribution n q (r) at r = 0.6. (b) Total number 
of atoms Ncr/l A (full line) as a function of time, for 8 — —16. 
The normalization with respect to I 3 makes this quantity in- 
dependent of the quantization volume. The dashed line is the 
respective bosonic result, while the straight dotted line is the 
total number of atoms N a (T)/l s ~ y / ]5|r/(47r) obtained using 
Fermi's golden rule calculation of the molecular spontaneous 
decay rate |2C| in the linear regime. 



For comparison, in the case of bosons the last terms in 
Eqs. (|2a|l and (j3a(l acquire positive signs, while the sin 
(cos) terms in the coefficients Aq and -B q are replaced by 
sinh (cosh) together with g q = (1 — (5 2 ) 1 / 2 and |A q | 2 - 
B q = 1 [23|. As a result, the solutions give n q (r) = 
sinh 2 ( 5q r) /.g 2 , and |m q (r)| 2 = [1 + n q (r)]n q (r). 

From Eq. (J3J we see that the momentum distribution 
of the atoms in the two spin states, f and |, is the same. 
The average mode occupancies undergo oscillations char- 
acteristic of fermionic statistics (see Fig. QJi); the maxi- 
mum occupancy of n q (r) = 1 imposed by the Pauli ex- 
clusion principle is reached at integer multiples of time 
t = 7r/2, for resonant modes satisfying g qo = 1. For qo to 
be nonzero, this condition requires a negative detuning 5, 
and therefore the absolute resonant momentum is given 
by go = |qo| = vRl- During the initial stage (r < 0.6), 
the occupancies grow in-phase and the 3D momentum 
distribution is peaked on the surface of a spherical shell 
of radius go as shown in the inset of Fig. At later 
times the oscillations dephase and the distribution func- 
tion becomes more complicated in structure. 

The total number of atoms in each spin state, N a (r) — 
Eq n q( T )i as a function of time is shown in Fig. \\]p. 
The initial growth of N„(t) saturates at r ~ 1.2, after 
which we see non-trivial oscillations. This is a combined 
effect of Pauli blocking and the oscillatory behavior of 
the individual mode occupancies. We emphasize that 
the saturation in this model is obtained within the unde- 
pleted molecular field approximation, and thus is purely 
a consequence of Fermi statistics. By comparison, the 
same approximation for bosons leads to an exponentially 
growing output due to bosonic stimulation and hence to 
unphysical results in the long time limit. Here, once the 
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FIG. 2: Snapshots of the average column density in momen- 
tum space n p (r)/l at different times r, for S = —16. 



depletion is taken into account, the saturation is natu- 
rally reached due to a finite initial number of molecules. 

In this respect, the present fermionic results suggest 
that the undcplctcd molecular field approximation for 
fermions is more reliable than for bosons. At the level 
of pairwise mode coupling, this conjecture is in fact sup- 
ported by the results of a numerical simulation using an 
exact quantum Monte Carlo method 28]. For example, 
at t ~ 7r/2 the discrepancy between the present and the 
exact results is about 5% for fermions, while it is ~ 11% 
for bosons and grows further with time. Similarly, the 
negligible role of the molecular depletion in the short time 
limit can be verified using the results of a related numeri- 
cal study of Ref. [2j| where the molecular field dynamics 
is treated at the mean-field level. 

We now turn to the analysis of pair correlations of the 
atoms in the opposite spin states and consider a normal- 
ized correlation function between 3D density fluctuations 
Ah n „ = n n 



= fiq a — (Hq CT ) in momentum space: 

<m(<W' T ) = ( Aft qT A "q'i>/\/ ("qtX'Vi! 



(7) 



For atom pairs with non-opposite momenta the pair 
correlation vanishes, 3fx(q, q', t)| ,^ = 0, implying 
the absence of any correlation. In the case of equal but 
opposite momenta, we find that 

ffufa -<L r) = |m q (r)| 2 /« q (r) = 1 - n q (r) < 1. (8) 

This corresponds to the maximum degree of correlation 
- as a consequence of Eq. ©, except when n q (r) = 1 in 
which case <?i-j.(q, — q, t) coincides with the uncorrelated 
level. For bosons, the respective pair correlation is given 
by </f j.(q, — q, r) = 1 + n q (r), which increases with n q (r) 
and always stays above zero. 

In order to make a better connection with the experi- 
ments at JILA Q] we note that the correlation measure- 
ments were made using absorption images after a time-of- 
flight expansion. This corresponds to analyzing the spa- 
tial column densities which involve integration of the 3D 
density along the direction of propagation of the imaging 
laser. Accordingly, we now analyze the momentum space 
analog of this procedure and calculate the correlation be- 
tween momentum column density fluctuations. 
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FIG. 3: (a) Correlation between momentum column density 
fluctuations, (p, — p,r), as a function of the absolute 2D 
momentum p — |p| at different times r, for 8 = —16. The 
correlation is larger initially, in the few-particle regime, and 
decreases as the number of atoms grows. Similar relationship 
is seen within the low-density tails and the higher-density 
central part of the momentum distribution, (b) Correlation 
signal as a function of time, for two different values of the p. 



The atom number operator corresponding to a z- 
integrated momentum column density is given by n pa — 
^qo- ; where p ={q x ,q y ) is the reduced 2D momen- 
tum. Using Eq. (J2J , the average column density is found 
via n p (r) = (n po -(r)) = ^ n q( T ) an d is the same for 
the two spin states. Snapshots of n p (r) at different times 
r are plotted in Fig. [5] The last two frames show a clear 
ring structure around the central background, which is 
consistent with the observed absorption images Q of spa- 
tial column densities after free expansion. 

The correlation function between momentum column 
density fluctuations in the two spin states, which we de- 
note via 5||(p, p', t), is defined as in Eq. (JTJ) except that 
the operators n qtT are replaced by n pa . In this case, the 
bars above the operators signify the procedure of summa- 
tion over the ^-component, before taking the ensemble 
average. For equal but opposite momenta, p' = — p, we 
find (see Fig. |3J) that 



5u(P>-P; r ) = 1 



E 



1z 



[nqMf/S P (r)<l. (9) 



For any other pair of momenta (p' ^ — p), the correlation 
function is simply zero, implying the absence of any cor- 
relation. In the case of dissociation into bosonic atoms, 
the same correlation function is given by g||(p, — p, r) = 

1 + E 9f K(r)} 2 /n p (r) >1. 

It is important to point out that the degree of corre- 
lation between atom pairs with opposite spins and mo- 
menta obtained in this model is maximal at any given 
density. The notion of maximal is defined here to corre- 
spond to perfect (100%) noise reduction of the number- 
difference fluctuations below the shot-noise level. 

This can be easily understood at the level of just two 
modes c q | and c_ q |, which we define via C\ and cV Con- 
sidering the normalized variance of the particle number- 
difference fluctuations, V = ([A(ni — fi 2 )\ 2 ) /SN, where 
SN is the shot noise level, one can show that in the sim- 
plest case of (ni) = (^2), the variance V and the pair cor- 
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relation g±2 = (Ani An 2 ) / \J (hi) (h 2 ) are related by V = 
l — 2gi2 (ni) / 5 AT. Here the shot noise level for fermions is 
given by SN = £\ ((An,) 2 ) = £. (h % ) (1 - (hi)), which 
we point out is always sub-Poissonian and is independent 
of the state of the fermion system. For a nonzero SN, the 
assumption of a perfect noise reduction below the shot- 
noise level, V = 0, can be used to identify the maximum 
degree of correlation, giving g^ ax ^ = 1 — (hi). This is 
exactly as obtained in the present model, Eq. (JBJ). At 
times when the average occupancies approach one, the 
shot-noise itself vanishes. Therefore, the notions of sub 
shot-noise fluctuations and maximal correlation become 
meaningless for (hi) = 1. 

For comparison, in the bosonic case the shot-noise level 
corresponds to that of a coherent state, SN = (hi) + 
(h 2 ), and never vanishes. For equal mode occupancies, 
it gives V = 1 + (hi) — 512 and therefore V — implies 
that g^ a ^ = 1 + (hi). This is always larger than the 
uncorrelated level of 0, and again agrees with the actual 
solution to the problem of molecule dissociation. 

In order to apply the results of the present uniform 
model to realistic trapped condensates, we remark that 
the quantization length L should be matched to the char- 
acteristic size of the molecular BEC. In addition, one has 
to ensure that the time window for dissociation is chosen 
such that the momentum kick fco — y/2m\A\/H imparted 
on the atoms is not too large, so that the atoms created 
mostly near the trap center remain within the molecular 
BEC while the dissociation is on. This implies that the 
present results are applicable to t < t max — Lm/(2hko) 
or T r max = l/(4y/\5\). Considering a typical set of 
paramet ers, with 3 x 10 5 initial number of molecules and 
I ~ 46 [23, the example presented here for 5 = —16 
would produce ~ 1.5 x 10 4 atoms in each spin state at 
t = 0.6 which compares favorably with r max ~ 3. This 
corresponds to 5% conversion and is consistent with the 
use of the undepleted molecular field approximation. 

In practice, the main factors which may contribute to 
the reduction of the correlation signal are: (?) the pres- 
ence of a large thermal component in the initial molecular 
gas, in which case the thermal centre-of-mass momenta 
may no longer be negligible compared to the momentum 
kick of the atoms ko; and (ii) atom-atom s-wave scatter- 
ing between the two spin components, which at high den- 
sities may substantially redistribute the momenta over 
the s-wave scattering spheres and spoil the correlations. 

In summary, we have analyzed short time dynamics 
of dissociation of a BEC of molecular dimers into corre- 
lated fcrmionic atoms in two different spin states. The 
pair correlations between atoms with opposite spins and 
momenta calculated here correspond to the maximum 
possible degree of correlation and serve as upper bounds 
for more detailed calculations and comparisons with ex- 
periments. The system may find applications in precision 
measurements beyond the shot-noise level, as well as for 



fundamental tests of quantum mechanics with macro- 
scopic number of fermions, such as demonstrations of 
Bohm's version of the Einstein-Podolsky-Rosen paradox 
[Tij and tests of Bell's inequalities for spin observables. 
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